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Abstract 



We have searched for the charmless hadronic decay of B° mesons into two 
neutral pions. Using 9.13 fb _1 taken at the T(4S) with the CLEO detector, we 
obtain an improved upper limit for the branching fraction B(B° — ► 7r°7r°) < 
5.7 x 10~ 6 at the 90% confidence level. 
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CP violation in the neutral kaon system was observed long ago JTJ, and evidence of 
CP violation in B decays is beginning to be seen ||. In the standard model, CP violation 
arises naturally from a single complex phase in the Cabibbo-Kobayashi-Maskawa (CKM) 
quark- mixing matrix ||. Observation of the time-dependent CP- violating asymmetry in the 
decay B° — > 7r + 7r~ (charge-conjugate modes are implied throughout this Letter) would, in 
principle, give a measurement of the sum of the CKM phases (3 = arg(V^) and 7 = axg(V* b ). 
However, difficulties arise from the fact that the tree and penguin contributions to the 
B° — > 7r + 7r~ decay enter with similar amplitude and unknown relative phase. It is known 
already from the large ratio of branching fractions B(B — > K + n~)/B(B — > 7r + 7r~) that the 
penguin contribution is large The tree and penguin contributions can (in principle) be 
separated by performing an isospin analysis on the related B — ► ttti modes ||. Although 
the decay mode B° — > 7r + 7r~ has been observed and there is some indication for the mode 
B + — > 7T + 7r° ||, the B° — > 7r°7r° decay mode has not been seen yet. Theoretical calculations 
offer possible values for the branching fraction B(B° — > 7r°7r°) ranging from 10 -8 to 10~ 5 JT]. 

In this Letter, we present a new limit on the B° — > 7r°7r° branching fraction based on data 
taken with the CLEO II detector at the Cornell Electron Storage Ring (CESR). The data 
consist of 9.13 ffcT 1 taken at the T(4S), corresponding to 9.67 x 10 6 BB pairs,[] and 4.35 ftT 1 
taken below the BB threshold, used for background studies. The new result supersedes the 
result from a previous publication ||, which was obtained with one third of the present 
statistics. 

CLEO II is a general purpose detector, described in detail elsewhere 0. Most relevant for 
the present analysis are the tracking system and the electromagnetic calorimeter. Momentum 
and specific ionization (dE/dx) of charged tracks are measured in cylindrical drift chambers 
in a 1.5 T solenoidal magnetic field. In a second configuration of the detector, CLEO II. V, the 
innermost tracking chamber was replaced by a 3-layer, double-sided silicon vertex detector, 
and the gas in the main drift chamber was changed from argon-ethane to a helium-propane 
mixture. These modifications improved both the charged particle momentum resolution 
and the dE/dx resolution. Photons are detected using a high resolution crystal CsI(Tl) 
electromagnetic calorimeter, composed of 7800-CsI(Tl) crystals covering 96% of the solid 
angle. Approximately two thirds of the data used in the present analysis was taken with the 
CLEO II. V detector. 

Photon candidates are defined as isolated showers with energy greater than 30 MeV in 
the central region of the calorimeter (| cos^l < 0.71 , where 9 is the polar angle relative to the 
beam axis), and with energy greater than 50 MeV elsewhere. Neutral pions are formed from 
pairs of isolated photons with invariant mass within 2.5 a of the 7r° mass (a rs 9MeV/c 2 
for a 2.7GeV/c 7r°). We require at least one of the two photons forming a 7r° candidate to 
be in the central region of the calorimeter. The energies of the selected photons are then 
kinematically fitted with the mass constrained to the 7r° mass. 

The B decay candidates are selected using a beam-constrained B mass M = \J E% — p B , 
where -Ebeam is the beam energy and pb is the B candidate momentum, and an energy 
difference AE = E\ + E 2 — -Ebeam, where E% and E 2 are the energies of the two neutral 



1 We assume equal branching fractions for T(45) — > B°B° and B + B 
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pions. The resolution in M is about 3.4MeV/c 2 , due to equal contributions from the beam 
energy spread and the 7r° energy resolution. The resolution on AE is approximately 60 MeV 
and is slightly asymmetric because of energy loss out of the back of the Csl crystals. Using 
events containing at least three charged tracks, we select B candidates with M in the range 
5.2 - 5.3 GeV/c 2 and \AE\ < 400 MeV. The fiducial region in M and AE includes the signal 
region and a substantial sideband used for background normalization. 

The main background arises from e + e~ — > qq (where q = u,d, s, c). Such events typically 
exhibit a two-jet structure and can produce high momentum back-to-back particles (tracks 
and/or showers) in the fiducial region. To reduce contamination from these events, we 
calculate the angle 6s between the sphericity axis |[D] of the candidate showers and the 
sphericity axis of the rest of the event. The distribution of cos 6s is strongly peaked at ±1 
for qq events, due to their two-jet structure, and is nearly flat for BB events. We require 
| cos^l < 0.8 which eliminates 85% of the qq background. Using a detailed GEANT-based 
Monte Carlo simulation JTTJ we determined the overall B° — > 7r°7r° signal detection efficiency 
£ = 28.8 %, dominated by geometric acceptance. Additional discrimination between signal 
and qq background (continuum) is obtained from event shape information used in a Fisher 
discriminant (J 7 ) technique as described in detail in Ref. [12| . 
used in the maximum-likelihood fit described below. 



The Fisher discriminant is 



A total of 1134 £ c 



7T°7r 



candidates are selected with the requirements described above. 
Figure [I] shows the AE versus M distribution of all the selected events, and the individual 
distributions of M, AE and with restrictions on the other two variables to emphasize the 
signal region. 
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FIG. 1. Distributions of the selected events, a) AE versus M for all selected events (only 
|cos#s| < 0.8 restriction is applied). The solid lines show the 2a boundaries, b) M distribution 
after 2a requirements on AE and T < 0.6. c) AE distribution after 2a requirements on M and 
T < 0.6. d) T after 2a requirements on AE and M. The dotted line shows the position of the cut 
on T for the other two variables. In plots b), c) and d), the solid line shows the result of the fit for 
signal and background, and the dashed line shows the contribution of the background alone. 



Monte Carlo simulation was also used to study backgrounds from b —* c and other b 
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and b — > s decays. More than 40 decay modes of the B meson into final states containing 
energetic tt°s and/or photons were considered. Only the B + — > p + 7r° decay channel was 
found to give a non-negligible contribution to the selection of B° — > 7r°7r° signal events. 
The three-body final state 7r + 7r°7r° can be misidentified as a two-body -k°7t° signal candidate 
when the charged pion from the asymmetric decay of the polarized p + has sufficiently low 
momentum and the ir° energies are poorly measured. The best separation between this 
background and signal is obtained in the AE distribution. The B + — > p + 7r° background is 
accounted for in the maximum-likelihood fit as described below. 

To extract the signal yield, we perform an unbinned maximum-likelihood fit using the 
variables M, AE, and T for each candidate event. The likelihood of an event is parameterized 
by the sum of probabilities of signal, qq background, and B + — > p + 7r° background hypotheses, 
with relative weights determined by maximizing the likelihood function C. The probability 
of a particular hypothesis is calculated as the product of the probability density functions 
(PDFs) for each of the input variables. Further details about the likelihood fit can be found 
in Ref. [fL2| . The PDFs for signal and B + — > p + 7i° are determined from high-statistics 
Monte Carlo samples. The PDFs for continuum are obtained from the data taken below the 
BB threshold. 

Monte Carlo experiments are generated to test the fitting procedure, and to produce 
frequentist confidence intervals as defined in Ref. |1|]. We generate Monte Carlo sam- 
ples containing the same number of events as the real data sample. Continuum events 
are generated according to the continuum PDFs, neglecting the small correlation between 
the fit variables. According to a given branching fraction B(B° — > 7r°7r°), we include sig- 
nal events randomly selected from our large B — > 7r°7r° Monte Carlo simulated sample. 
We also include Monte Carlo simulated B + — > p + 7r° events. We generate 1000 samples 
for several values of B(B° — > 7r°7r°) and B(B + — > p + 7r°) in the range — 10 x 10~ 6 
and — 42 x 10~ 6 , respectively. We apply the fitting procedure to every sample indi- 
vidually and determine the signal yield distribution for each value of B(B° — > 7r°7r°) and 
B{B + — > p + 7T°). In the samples containing B + — > p + 7r° events, we find a small increase of 
the signal yield proportional to B(B + — ► p + 7r°). At the 90% confidence level (C.L.) upper 
limit B(B + —> p + n°) < 42 x 10~ 6 flI5] , the contribution to the signal yield is 0.3 event. We 
include this maximal contribution as a one-sided systematic uncertainty in the result. The 
Monte Carlo experiments show that once the B + — ► p + 7r° background is accounted for, the 
average yield for any value of B(B° — ► 7r°7r°) is equal to the expected yield for this branching 
fraction, excluding thus any significant bias from the fitting method. An estimation of the 
statistical sensitivity of the measurement is given by the width of the yield distributions, 
and is measured to be about ±5 events. 

Figure shows the 68%, 90%, 95%, and 99% frequentist confidence intervals (statistical 
only) built from the signal yield distributions obtained with B{B + — > p + 7r°) = and following 
the method described in Ref. ||14|| . 

Figure |3] shows the result of the likelihood fit as a plot of x 2 ~ Xmm — ~ 2 ln£/£ m ax- ^ ne 
maximum likelihood £ max is found for a signal yield N$ = 6.2l 3 ; 7 events, with a statistical 
significance of 2.0 a. We define the statistical significance to be no if the value of — 21n£ 
increases by n 2 when the signal yield N$ is constrained to be zero. The measured yield for 
B + — > p + 7r° is N pn = — 11 ± 9 events, consistent with the upper limit for that mode [|i~3| . 



This yield should, however, not be used to calculate a new value of the upper limit for 



6 




) versus signal 



yield Ns, as determined from Monte Carlo experiments. No systematic effects are included. The 
vertical line indicates the likelihood fit result for Ns- 

B{B + — > p + 7r°), as the 7r°7r° analysis described here is deliberately designed to minimize 
sensitivity to the B + — > p + 7r° mode. We also checked that, due to the small correlation 
between signal and background, the signal yield is reduced by only 0.2 event when the 
background yield is constrained to be positive. Given the measured value of Ns, we use 
Figure |2] to determine the frequentist confidence interval for B(B° — > 7r°7r°). We obtain the 
90% C.L. interval 0.3 x f0~ 6 < B(B° -> 7r°7r°) < 5.2 x 10" 6 . This interval is statistical only, 
and does not include systematic uncertainties. 

For the treatment of systematic uncertainties, we separate them into two categories. 
First, we estimate a systematic uncertainty on the fitted signal yield by varying the PDFs 
used in the fit within their uncertainties, and we add in quadrature the previously described 
systematic uncertainty due to the possible residual contamination from B + — ► p + 7r° events. 
We also consider a possible mismodeling of the e + e~ — > t + t~ contribution in the PDFs for 
continuum, and we assign a systematic uncertainty of ±0.7 event based on a high-statistics 
Monte Carlo simulation study. We obtain the total systematic uncertainty t{_ 8 events on the 
signal yield Ns- Secondly, we estimate an uncertainty on the signal detection efficiency £ , to 
account for uncertainties related to 7r° finding efficiency, maximum-likelihood fit efficiency, 
luminosity, Monte Carlo statistics, and the |cos#s| requirement. The efficiency with its 
uncertainty is £ = (28.8 ± 2.0)%. 

We derive the central value of the branching fraction B(B° — > 7r 7r°) = ^2.2 jig £J x 
f 0~ 6 , where the first uncertainty is statistical and the second is systematic. We also calculate 
the 90% C.L. upper limit yield by integrating the likelihood function 
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g C m UN)dN 
I °° jC max (N) dN U ' yU ' 

where £ max is the maximum £ at fixed N to conservatively account for possible correlations 
among the free parameters in the fit. We obtain the upper limit N s < 13.7 events at 90% 
C.L. (statistical). We then calculate the corresponding upper limit of the branching fraction, 
add one standard deviation of the systematic uncertainty, and obtain the branching fraction 
90% C.L. upper limit B(B° -> 7r°7r°) < 5.7 x 10~ 6 . 

In summary, using the full CLEO II and CLEO II. V data set, we have obtained an 
improved upper limit on the branching fraction of the B° — > 7r°7r° decay mode. We see no 
indication for a signal and we set a new 90% C.L. upper limit B(B° — > 7r°7r°) < 5.7 x 10~ 6 . 
This limit is in the range of the theoretical predictions and constrains some of them. 
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